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ABSTRACT 



Two basic formation channels have been proposed for blue straggler stars in glob- 
O I ular clusters: binary star evolution and stellar collisions. We recently showed that the 

number of blue stragglers found in the core of a globular cluster is strongly corre- 
^ ■ lated with the total stellar mass of the core, but not with the collision rate in the 

core. This result strongly favoured binary evolution as the dominant channel for blue 
straggler formation. Here, we use newly available empirical binary fractions for glob- 
ular clusters to carry out a more direct test of the binary evolution hypothesis, but 
^ . also of coUisional channels that involve binary stars. More specifically, using the cor- 

relation between blue straggler numbers and core mass as a benchmark, we test for 
correlations with the number of binary stars, as well as with the rates of single-single, 
single-binary, and binary-binary encounters. We also consider joint models, in which 
blue straggler numbers are allowed to depend on star/binary numbers and collision 
I rates simultaneously. 

Surprisingly, we find that the simple correlation with core mass remains by far the 
strongest predictor of blue straggler population size, even in our joint models. This is 
despite the fact that the binary fractions themselves strongly anti-correlate with core 
mass, just as expected in the binary evolution model. 
_ At first sight, these results do not fit neatly with either binary evolution or col- 

rS , lisional models in their simplest forms. Arguably the simplest and most intriguing 

■ possibility to explain this unexpected result is that observational errors on the core 

binary fractions are larger than the true intrinsic dispersion associated with their 
dependence on core mass. In the context of the binary evolution model, this would 
explain why the combination of binary fraction and core mass is a poorer predictor of 
blue straggler numbers than core mass alone. It would also imply that core mass is a 
remarkably clean predictor of core binary fractions. This would be of considerable im- 
portance for the dynamical evolution of globular clusters, and provides an important 
benchmark for models attempting to understand their present-day properties. 

Key words: blue stragglers - binaries: close - globular clusters: general - methods: 
statistical - stellar dynamics. 



1 INTRODUCTION 
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Blue stragglers (BSs) in globular clusters (GCs) are stars 
that appear brighter and bluer than the main-sequence 
turn-off (MSTO) in the cluster colour-magnitude diagram 
(CMD) (jSandagd[l953 ). They are thought to be created 
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when fresh hydrogen is mixed into the core of a nor- 
mal l ow-m ass main-sequence (MS) star (e.g. ISills et al] 
I2OOII . I2OO2I '). Two main formation channels have been pro- 
posed for BSs in GCs: binary star evolution and dynami- 
cal interactions. BSs can form via the former pathway ei- 
ther through binary mass-transfer due to Roche-lobe over- 
flow fro m an evolved primary onto a norm al MS com- 
panion (|McCrealll963 : lGeller fc Mathieull201ll '). or through 
the coalescence of two normal MS stars in a binary sys- 
tem. The latter mechanism can occur due to angular mo- 
mentum loss induced by a magnetize d stellar w ind (e.g. 
Iben fc Tutukovlll984l : lAndronov, Pinsonneault fc Terndrud 



20061 ), or e ven Kozai cycles induc ed by an outer triple 



companion (|Perets fc FabrvckvllioogI ') . The dynamical path- 
way involves collisio ns between two or more MS stars (e.g. 
ISills fc Bail^ll999l ). These are typically mediated by dy- 
namical interactions involving binary stars, since the cross- 
section for collisi on is much larger for a binary than it is for 
a smg le star fe.g. iLeonardll 19891 : [Leonard fc Linnelllll992l ). 

Several statistical studies have been conducted in 
search of a dominant BS formation channel. However, 
the cluster parameter that has thus far yielded the 
stro ngest correlation with BS population size is t he clus- 
' Knigge '2007*) or core 



ter (iPiotto et al. 2004 



Lei gh. Sills 
200' 



l|Knigge. Leigh fc S ills 2009|; iLeigh. Sills fc Kniggj l2011al ') 
mass. Many authors have tried to explain this by simultane- 
ousl y invoking multiple formation mec hanisms. For exam- 
ple, bavies. Piotto fc De Angelil (|2004 ) suggested that the 
observed dependence of BS numbers on total cluster mass 
can be explained if BSs in low-mass clusters are primarily 
descended from binaries, whereas BSs in high-mass clusters 
are primarily descended from collisions. A similar scenario 
has been argued for in an attempt to explain the bimodal BS 
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BSs in the dense core were formed in collisions, whereas BSs 
in the low-density cluster outskirts were formed by mass- 
transfer within primordial binaries. 

At the time these studies were conducted, there were 
hardly any observational constraints on the properties of the 
binary populations in GCs. In particular, empirical binary 
fracti ons were available fo r only a small subset of low-density 
GCs llSollima et al. 200S ). Th is issue was only recently re- 
solved bv lMilone et al.l l|2012h . Using data from the HST- 
based ACS Survey of Globular Clusters, these authors de- 
rived photometric binary fractions for the MS populations in 
59 Milky Way (MW) GCs. This sample offers a long-awaited 
opportunity to test more directly whether the sizes of BS and 
binary populations in GCs are correlated, as one might ex- 
pect for both the binary evolution channel and for coUisional 
formation channels that involve binaries (i.e. l-|-2 and 2-|-2 
encou nters) . This was previousl y addressed bv lSoUima et al.l 
l|2008h and iMilone et al.1 (|2012l '). These results provided ev- 
idence that the blue straggler fraction is indeed related to 
the binary fraction, and we build on those previous works 
in this paper. 

The plan of this paper is as follows. First, in Section [2l 
we present the observational data. We then explain and 
carry out our analysis of these data in Section |31 There, we 
derive the expected scaling laws for the simplest versions of 
the various formation channels and compare these theoreti- 



cal predictions to the observations. Finally, the implications 
of our results for the formation and evolution of BSs in GCs 
are discussed in Section 3] 



2 OBSERVATIONAL DATA 

BS num bers are taken from Table 1 of lLeigh. Sills fc Kniggel 
(|2011ah . which was compiled u sing data taken from th e ACS 
Survey for Globular Clusters jSaraiedini et al.l l2007l f1 The 
sample used in this paper omits five clusters from the cata- 
logue of Leigh, Sills fc Knigge (20 11ail , since we do not have 
observed binary fractions in these cases. We use only those 
BS number counts within the core and within four core radii 
from th e cluster centre (columns 4 and 7, respectively, in Ta- 
ble 1 of lLeigh. Sills fc Knigg3 l|2011ah ') in this paper. 

Binary fractions within the core (r < rc5 ffein,core) and 
the half-mass radi u s (r < r^; ibin,haif) are taken from Table 
1 of iMilone et al l (I2OI2I). B T h e latt er values are not pro- 
vided directlv in iMilone et al.l l|2012t ). Instead, binary frac- 
tions within the annulus separating the core and the half- 
mass radius (rc < r < r^; ibin,r^<r<r^) are given. Therefore, 
we calculate mass-weighted binary fractions within the half- 
mass radius using the relation: 

[Mhalf - Mcore) 



M, 



fbi n,half 

(1) 

where Mcore and M^haif are the mass of the cluster core 
and the mass contained within the half-mass radius, re- 
spectively. In order to obtain accurate estimates for the 
total stellar mass contained within the core, we gener- 
ated sin gle-mass King models c alculated using the method 
of .Sigurdsson fc PhinnevI (|l995l ) to obtain luminosity den- 
sity profiles for every cluster in our sample. The pro- 
files were found using the c oncen tration parameters of 
iMcLaughlin fc van den Marell (|2005|) and the central lumi- 
nositv densities of iHarriX i 1996. 2010 updatel) . We then in- 
tegrated the derived luminosity density profiles numerically 
in order to estimate the total stellar light contained within 
the core, which we multiplied by a mass-to-light ratio of 2 
in order to obtain estimates for the total stellar mass con- 
tained within one core radius from the cluster centre. The 
mass enclosed within the half-mass radius was then esti- 
mated by calculating the total cluster mass from its absolute 
integrated visual magnitude (once again assuming a mass- 
to-light ratio of 2), and then dividing by two. Throughout 
this paper, we have ad opted the bina r y frac tions provided in 
column 6 of Table 1 in IMilone et all |2012f ). which provides 
the total fraction of objects that are binaries within the in- 
dicated annulus. However, in practice, this estimate of the 
total number of binaries is extrapolated from the observed 
fraction of binaries with mass ratio q > 0.5 by assuming a 
flat distribution in q. 



^ The data can be found 

http: / /www. astro. ufl. edu/^ata/public-hstgc/[ and 



at 

„ , was last 

accessed on 02/02/11. ' 
^ These binary fractions have been corrected for a variety of ob- 
servational biases, including completeness, contamination from 
field stars, and differential reddening. Detailed explanations o f 
these procedure s have bee n provided in ISaraiedini et all l[2003), 
lAnderson et al.l 1I2OO8I ') and IMilone et al.l ll20ia V 
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Figure 1. The binary fraction in the core plotted versus the bi- 
nary fraction in the annulus separating the core and the half-mass 
radius . Error bars are taken directly from Table 1 of Mi lone et al.l 
(120121') . The solid line shows the weighted least-squares fit to the 
data provided in Equation [2] 



The re are severa l clust ers for which the binary frac- 
tions in iMilone et al.l ()2012l ) are provided for f(,in,rc<r<rfc 
but not ibin,core- To approximate fbin,core in those clusters 
for which these values are missing, we perform a weighted 
least-squares fit to quantify the dependence of ibm.corE on 
fbin.r, <r<r^ using evcry cluster in the sample of ! Milone et al.l 
l|2012l) for which both of these quantities were given. We 
then supplement ibin,core in every cluster for which only 
ibin,rc<r<r,^ was provided. This is only necessary in a hand- 
ful of clusters, but the resulting increase in our sample size is 
nevertheless worth while. We obtain a relation of the form: 

fb^n.core = (0.84 ± 0.15) fb^„,r, <r<r^ + (0.04 ± 0.01). (2) 

This is shown in Figure [1] The fit is good for binary fractions 
less than ~ 0.2, but the agreement is poor for larger binary 
fractions. In order to test the effects had on our least-squares 
fit by clusters with large binary fractions, we re-perform it 
considering only those clusters for which both fbin,core < 0.2 
and ibin,rc<r<r,^ < 0.2. In this case, both the slope and y- 
intercept of the fit agree with those reported in Equation [2] 
to within one standard deviation. We conclude that the fit 
and the corresponding uncertainties provided in Equation [2] 
are reasonable. 

Our analysis in Section |3] also requires the absolute 
visual magnitude (My), core radius (tc), central luminos- 
ity density (po) and central velocity dispersions (ap) of 
each cluster. All of these quantities are taken from iHarrij 
(|1996. 2010 updatel) . except for the velocity dispersion of 10 
clusters for which these values are not provided. In these 
ca ses, we use the calculated velocity dispersions provided 
bv lWebbinkI (|l985l ). 



3 ANALYSIS AND RESULTS 

In this section, we derive theoretical scaling laws for the 
simplest versions of the various proposed blue straggler for- 
mation channels and compare these theoretical predictions 
to the observational data. In practice, this means we will 
search for correlations between the number of blue strag- 
glers in the core of each GC and the cluster parameter that 
should set this number in each scenario. 

We assess the significance of our correlations primarily 
via the Spearman rank test, which provides both a correla- 
tion coefficient (ts) and the significance level at which the 
null hypothesis of zero correlation is disproved (ps). A small 
Ps-value is indicative of a significant correlation. In addi- 
tion, we perform weighted least-squares fits to quantify the 
dependence of BS numbers on each parameter. Uncertainties 
for all number counts are obtained assuming Poisson statis- 
tics, but we also allow for intrinsic dispersion in our fits, at 
whatever level is required to achieve a reduced — 1- 

Table [T] shows the results of our comparisons between 
the observed BS numbers and the parameter we have tested. 
Each entry in Table [T] gives the slope for the line of best-fit. 
Spearman correlation coefficient, and probability at which 
the null hypothesis of zero correlation is disproved. These 
are provided in the form (slope; rs, Ps). 

3.1 Core Mass 

We begin by revisiting the correlation between core BS num- 
bers and core mass, with the latter estimated from our King 
model fits. Core mass was fo und to be the best pr e dictor 
of BS numbers in the core in iKnigge. Leigh fc Silld (|2009l ) 
and lLeigh fc Silli l|2011bl ). Our working assumption in those 
studies was that core mass was able to predict BS numbers 
because it is a partial proxy for the (unknown) number of 
binaries in the core (since Nbin,core = f bin, cor eM core', see 
Section r3.2|l . Now that binary fractions are available, we can 
test this assumption directly. However, the original correla- 
tion with just core mass remains a useful benchmark in this 
context: if the empirical binary fractions have added useful 
information, including them should allow us to discover even 
stronger correlations. 

We find a dependence Nss.core oc M^o^e^"'"^, and a 
Spearman correlation coefficient 0.83. This correlation is 
shown in Figured along with the corresponding line of best- 
fit to the data. The power-law index on Mcore is inconsis- 
tent with zero at the 8a confidence level, but also inconsis- 
tent with unity at the 12a confidence level. Therefore, the 
strong, sub-linea r correlation between NBS.co rc and Mcore 
first reported in iKnigge. Leigh fc Sillsl (|2009l ) is confirmed 
in the present analysis. 

3.2 Binary Population Size 

If most of the BSs in our sample are descended from binary 
evolution, we predict a dependence of the form: 

jvr TIT fbin,rMencL,r /q\ 
iyBS.core OC I\bin,r ~ ' , (o) 

m 

where Nbs is the number of BSs within a given radius r from 
the cluster centre, Nbin,r is the number of binaries contained 
within r, fbin,r is the fraction of objects within r that are 
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Figure 2. The logarithm of the number of BSs as a function of 
the core mass (in units of solar masses). The solid line shows the 
best-fit to the data. 



binaries, Menci,r is the total stellar mass enclosed within r, 
and fh is the average stellar mass (for which we assume the 

same value in all clust ers) . 

As pointed out in lKnigge. Leigh fc Sillj l|2009l ). Equa- 
tion |3] allows a simple explanation for the observed sub- 
linear correlation between Nss.core and Mcore- If we identify 
Mcore with Menci,r (i-S- wc take the relevant radius to be 
r = rcore), then Nss.core oc Mcore follows from Equation [3] 
if the core binary fraction is itself a function of core mass, i.e. 

i 'bin. co re OC M~Jlf . We also showed in , Knigge. Leigh fc Sillsl 
20091 ) that there was indeed some evidence for an anti- 
correlation between cor e mass and core bin ary frequency, 
based on the data set of ISollinia et al.l (|2008l ): this IS a pre- 
limimary version of the ACS-based data set we use here. 

The final ACS binary fraction data set available now 
iMilone et all I2012D allows us to check whether the evi- 
dence for this anti-correlation holds up. Figure [3] shows that 
it does. More specifically, we find fbin.core oc M^J^f^'^'''^, 
while the Spearman correlation coefficient is r^ = —0.72. 
The slope for this anti-correlation is not quite as steep as in 
our naive prediction, but the existence of the anti-correlation 
itself is obviously promising. 

However, this promise is not actually borne out. When 

we compare NnS.core directly to NMn,core = fbin,coreMcore 

(Figure [31 Table [1]), we find that the strength of the corre- 
lation actually decreases compared to the correlation with 
just Mcore (the Spearmau rank coefficient drops from 0.83 
to 0.61), while the best-fit slope increases only marginally 
(from 0.40 ± 0.05 to 0.48 ± 0.09). Increasing the size of the 
region considered does not improve things: when we compare 
the number of BSs within four core radii from the cluster 
centre against the number of binaries within the half-mass 
radius, we still find only a weak correlation with a clearly 
sub- linear slope (Figure [5] Table [T]). This is a surprising re- 
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Figure 3. The logarithm of the core binary fraction as a function 
of the core mass (in units of solar masses). The solid line shows 
the best-fit to the data. The seven data points with open circles 
around them correspond to the core binary fractions calculated 
using Equation [2] 

suit at first sight, and we will discuss its implications further 
in Section 131 

Before moving on, it is worth stressing that the strong 
anti-correlation between binary fraction and core masse 
shown in Figure|3]is interesting in its own right. Empirically, 
its existence is in line with a similarly strong correlation be- 
tween core binary fraction an d cluster absolute m agnitude 
that was already presented in lMilone et al. I l|2012 h). This is 
because absolute magnitude is a proxy for total cluster mass, 
which in turn correlates strongly with core mass among GCs. 
Theoretically, however, it is far from clear why abundance 
of binary stars should depend so strongly on either the core 
or the total stellar mass of their host clusters. The fact that 
it does must be important for our understanding of cluster 
dynamics. Presumably, dynamical evolution is responsible 
for establishing this correlation and is, in turn, affected by 
it. 

3.3 Collision Rates 

If most of the BSs in our sample were formed from direct 
collisions between single stars, then we predict a dependence 
of the form: 

Nbs cc Ni+i ~ / Fi+idt, (4) 

J TO 

where Fi+i — 1/ri+i is the rate of single- si ngle collisions 
produc ing BSs (we use the form given in iLeigh fc SilQ 
(|2011bl )). and we are integrating with respect to time. We 
use the age of the cluster (tci ) as the upper limit of integra- 
tion, and the average BS age {tbs) to calculate the lower 
limit of integration according to ro = Tci —tbs- Assum- 
ing Fi+i remains constant in time l|Leigh. Sills fc Kniggel 
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Figure 4. The logarithm of the number of BSs (within the core 
only) as a function of the number of binary stars in the core. The 
solid line shows the best-fit to the data (the slope for which is 
shown in Table [TJ- 

l201ld ). Equation |4] simplifies to: 

Nbs oc Ni+i ~ Ti+iTBS- (5) 

If most of the BSs in our sample were formed from dy- 
namical interactions involving binaries, then we predict ei- 
ther a relation of the form: 

Nbs oc A''i+2 ~ Ti+2rBS, (6) 

where ri-1-2 = l/ri+2 is the rate of single-binary encounters, 
or: 

Nbs oc A^2-i-2 ~ ^2+2rBS, (7) 

where r2-i-2 = I/T2-1-2 is the rate of binary- binary encounters. 

We find a correlation between the number of BSs in the 
core and the 1-1-1 collision rate. The slope for this relation 
is inconsistent with zero at nearly the 4cr confidence level, 
and it yields a Spearman correlation coefficient that is only 
slightly smaller than we find for the numbers of binaries in 
the core. For the 1+2 and especially the 2-1-2 collision rates, 
however, our results are consistent with little to no correla- 
tion with BS population size. The relations for the 1-1-1, 1+2, 
and 2-1-2 collision rates are plotted in Figure[6l FigureO and 
Figure [51 respectively, along with the corresponding lines of 
best-fit given in Table [T] 

3.4 Joint Models 

Several authors have suggested that multiple mechanisms 
could contribute significantly to BS formation. In particu- 
lar, different formation mechanisms c ould operate simulta - 
neously within the same cluster (e.g. iFerraro et al.' 2004h. 
and /or different format ion mechanisms could dominate i n 
different clusters (e.g. iDavies. Piotto fc De Angelil 12004 ). 




2.5 3 3.5 

log \" . , 



Figure 5. The logarithm of the number of BSs (within four core 
radii from the cluster centre) as a function of the number of binary 
stars within the half-mass radius. The solid line shows the best-fit 
to the data. 



One prediction that is consistent with the second scenario is 
that the BSs in the high-density clusters in our sample were 
formed from collisions, whereas the BSs in the low-density 
clusters were formed from binary star evolution. The first 
scenario, on the other hand, predicts that some linear combi- 
nation of the parameters Ni,i„, Ni+i, N1+2, and N2-1-2 should 
yield the strongest correlation with observed BS numbers. 

We perform two additional comparisons in an effort 
to test the idea that multiple formation mechanisms con- 
tributed to the formation of the BSs in our sample. First, we 
divide our sample into low- (logpo < 3.3) and high-density 
(log po > 3.3) sub-samples of roughly equal size, and inde- 
pendently re-perform our analysis on each. This comparison 
will tell us whether different formation mechanisms dom- 
inate in each of these sub-samples independently. We use 
the cluster density to divide our sample since, if collisions 
do contribute to BS formation, they should occur with the 
greatest frequency in high density clusters. 

For both sub-samples of low- and high-density clusters, 
our results remain consistent with what we found for the 
entire sample to within the uncertainties. This is supported 
both by our lines of best-fit and our Spearman correlation 
coefficients. Therefore, this is consistent with the general pic- 
ture that the dominant BS formation mechanism is the same 
in all clusters (as opposed to different mechanisms dominat- 
ing in different clusters). 

Next, we search for a linear combination of the different 
formation channels that yields a better correlation with the 
observed BS numbers than any of these parameters individ- 
ually. This is done in two ways. Specifically, we fit to the 
observed data relations of the form: 

Nbs = aNun,core + bNl+l + cNl + 2 + dN2 + 2, (8) 
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Figure 6. The logarithm of the number of BSs as a function of 
the single-single collision rate in the core (in units of number of 
collisions per year). The solid lines shows the best-fit to the data. 
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Figure 8. The logarithm of the number of BSs as a function of 
the binary-binary collision rate in the core (in units of number of 
collisions per year). The solid line shows the best-fit to the data. 
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Figure 7. The logarithm of the number of BSs as a function of 
the single-binary collision rate in the core (in units of number of 
collisions per year). The solid line shows the best-fit to the data. 

and 

NbS = eff^n,,oreM^ore + h fl^^^^rM+l . (9) 

where e, /, g, h, i, and j are all treated as free parameters, 
and we omit the factor (1 — fbin,core)~^ when calculating 
Ni+i from Equation [5] These comparisons will help to tell 



Table 1. Slopes for all weighted least-squares fits, and the corre- 
sponding Spearman rank correlation coefficients and probabilities 
that the null hypothesis of zero correlation is disproved. 



Parameter 


log Nbs 




slope; rs; ps 


log N&in,core 


0.48 ± 0.09; 0.61; 3.80e-4 


log ^bin,half 


0.53 ± 0.11; 0.49; 5.91e-3 


log Tl+i 


0.16 ± 0.04; 0.60; 5.14e-4 


log ri + 2 


0.13 ± 0.06; 0.36; 4.86e-2 


log r2+2 


0.04 ± 0.06; 0.06; 7.65e-l 


log Mcore 


0.40 ± 0.05; 0.83; 1.57e-8 



us if multiple formation mechanisms contribute significantly 
to BS formation in the same cluster, or if there is always a 
dominant formation channel. 

Equation [8] is unable to provide a more statistically 
significant correlation than we find between BS numbers 
and the core masses. The best-fitting model offers at best 
a slight improvement over what we find upon comparing 
BS population size to either the number of binaries in the 
core or the collision rate alone. Finally, the only best- 
fitting parameter in Equation |9] we find to be consistent 
with a non-zero value is the power-law index on core mass, 
with g = 0.48lo Jg. This independently confirms that the 
strongest dependence we find is between BS numbers and 
core mass. Therefore, we do not find evidence from this com- 
parison that multiple mechanisms contribute simultaneously 
to BS formation in individual clusters. 
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4 DISCUSSION 



4.1 Why do the empirical binary fractions fail to 
yield improved correlations for blue 
stragglers? 

Our main result is that core mass remains a better predictor 
of blue straggler numbers than any other variable we have 
considered. This is surprising. In the context of the binary 
evolution scenario for blue stragglers, one might have ex- 
pected a stronger correlation with the number of binaries 
in the core (as estimated by fbin,core x Mcore)- Similarly, 
in the context of the collision scenario, it seems reasonable 
to suppose that blue stragglers may form primarily in col- 
lisions involving binaries. In this case, stronger correlations 
with either 1+2 or 2-1-2 collision rates might be expected. 
Yet neither of these hypotheses is coniirmed by the data. 

At first sight, this negative result is all the more sur- 
prising because the binary fractions themselves do (anti- 
)correIate strongly with core m ass (and absolute cluster 
mass). As noted in lKnigge. Leigh fc Sil ls (200 9). this anti- 
correlation is exactly what is needed to explain the sub- 
linear correlation between BS numbers and core mass in the 
binary evolution model. So why does the combination of 
core mass and binary fractions not lead to a roughly linear 
correlation between binary and BS numbers, as one might 
naively expect? 

It is, of course, possible that the core mass correlation 
is simply the most fundamental one. However, if this cor- 
relation is not driven by binaries, its origin and sub-linear 
nature are quite hard to understand. 

As it turns out, there is, in fact, a simple way to un- 
derstand our results in the context of the binary evolution 
model. The easiest way to see this is to recognize that the 
existence of a correlation between M^ore and fbin,core (Fig- 
ure [H]) obviously implies that core mass is itself an estimator 
of core binary fractions. Thus even if Nbs depends solely 

on Nbin,core, replacing A/core with fbin,core X Mcore will Ouly 

lead to an improved correlation if the observational errors on 
the empirical binary fractions are smaller than the intrinsic 
scatter in the fbin,core versus Mcore relationship. 

In order to illustrate this quantitatively, we have carried 
out some simple simulations. Thus we create mock data sets 
containing = 30 data points and spanning roughly the 
same dynamical range as the real data. In each mock data 
set, we assume that the true number of BSs scales perfectly 
and linearly with the number of binaries. We also assume 
that binary fractions correlate sub- linearly with core masses, 
fbin,core oc M~Jlf , and that this correlation is quite tight, 
with an intrinsic dispersion of aint =0.1 dex. We also as- 
sume that BS numbers are subject to an observational error 
of 0.1 dex, and, for simplicity, that core masses are per- 
fectly known. Finally, we assume that the binary fractions 
are subject to observational uncertainties, aobs{fbin,core)- 
We are interested in how the character of observation- 
ally inferred correlations changes when (Jobs{fbin,core) ap- 
proaches and exceeds aint, so we run tests over the range 
Q.lGint ^ Oobsifbin.core) ^ IQ.Qoint- For each trial value of 
o'obs{fbin,core), wc Create 1000 mock data sets and measure 
the Spearman-rank correlation coefficients of the NBS,obs 
versus Mcore relation and the Nss.obs versus Nbi„,obs rela- 
tion. We also measure the slope of the log Nss.obs versus 
log Nbin,obs correlation in each mock data set, in order to 




S 0,2 - 



Figure 9. The results of our simulations to quantify wliethcr the 
observational errors on the empirical binary fractions are smaller 
than the intrinsic scatter in the core versus Mcore relation. 
Our procedure for this has been described in detail in the text. 
Note that the error bars shown in this plot correspond to the 
standard deviation across the mock samples, as opposed to the 
error for the mean. 



check if and when this deviates substantially from the true 
slope of unity. 

The results of the simulations are shown in Figure |9l 
They confirm that the correlation coefficient of the NBS,obs 
versus Nbin,obs relation exceeds that of the NBs,obs versus 



Mcore relation only if aobs {fk 



= ) 



Once the errors 



on fbin,core cxcecd this, the correlation with core mass be- 
comes stronger than that with Nbin,obs, as in the actual data. 
It is also interesting that, in the same regime, the measured 
slope of the \og Nss.obs versus log Nbin,obs relation becomes 
significantly shallower than the true slope of unity. Again, 
this matches what we see in our analysis of the actual data. 

We therefore suggest that the simplest way to under- 
stand our results is to assume that the observational un- 
certainties on the empirical binary fractions exceed the in- 
trinsic dispersion in the fbin,core versus Mcore relationship. 
If this is correct, the observational data can be understood 
in the context of the simple binary evolution model. It is 
worth stressing here that we are not suggesting that the 
observed binary fractions are "wrong" - merely that the un- 
certainties affecting them are larger than aint - This, in turn, 
is astrophysically important. The observed anti-correlation 
betwe en fbin.core and cor e mass (Figure [Sjl or total cluster 
mass (|Milone et al.l l|2012l) ) is already surprising and should 
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be a key benchmark for dynamical models of GCs. If our 
suggestion here is correct, this correlation is even tighter 
than the present data suggest. 

Having suggested our preferred explanation for our ini- 
tially surprising results, we will devote the rest of this section 
to explore the viability of other possibilities. In particular, 
we will carefully consider the potential impact of key as- 
sumptions in our analysis. 



4.2 A check on key assumptions in the analysis 

4-2.1 Constant average stellar mass 

We have assumed throughout our analysis that the average 
stellar mass, m, is constant across all clusters. In reality, 
the average main-sequence mass can range from ~ 0.3 M0 
to 0.7 Mq, and there seems to be a connection between 
mass function slope and the total GC mass, with lower-mass 
GCs being more depleted of preferentially low-mass stars 
(|Leigh et al.ll2012l l. This suggests that the average stellar 
mass should increase with decreasing cluster mass. How- 
ever, replacing the assumption of a constant average mass in 
Equation [3] with a dependence of the form m oc M^^^^, with 
e < 0, should further flatten the dependence of BS numbers 
on the calculated numbers of binaries. This suggests that 
our assumption for the average stellar mass is not the cause 
of the observed non-linear dependence of BS numbers on the 
numbers of binaries. 

We adopt the same assumption of a constant m when 
calculating the 1-1-1, 1-1-2, and 2-1-2 collision rates. The 
power-law indices we find with BS numbers for all three 
of these parameters are very small (~ 0.1), and the range in 
average stellar masses for our sample is at most ~ 0.4 Mq. 
Therefore, there is no realistic assumption for the average 
stellar mass that we could adopt to recover a linear relation 
between BS population size and any of the collision rates. 



4-2.2 Constant average semi-major axis for binaries 

Similarly, we assumed a constant value of a = 2 AU for the 
average semi-major axes of all binaries undergoing 1-f 2 and 
2-1-2 encounters. However, it is possible that d depends sys- 
tematically on the cluster mass, since the semi-major axis 
corresponding to the hard-soft boundary depends on the 
cluster mass (via the velocity dispersion). We replaced d — 2 
AU in our estimates for the 1-1-2 and 2-1-2 collision rates with 
the semi-major axis corresponding to the hard-soft bound- 
ary in each cluster, and re-performed our analysis for these 
two parameters. This did not improve the agrement between 
the observed BS numbers and either the 1-1-2 or 2-1-2 collision 
rates. 



4-2-3 Other cluster-to- cluster variations 

There are a number of ways that cluster-to-cluster variations 
in the distributions of binary orbital parameters could have 
affected our results. For instance, there are several reasons 
why high-mass MS-MS binaries with mass ratios g ~ 1 are 
the most likely to produce BSs. It is MS stars with masses 
just below that of the turn-off that are next in line to as- 
cend the giant branch. Provided they are in binaries, they 
are therefore the best candidates to over-fill their Roche 



lobes within the next few hundred Myrs. Two MS stars with 
masses close to the turn-off should als o produce the br ight- 
est and bluest BSs upon merging (e.g. ISills et al.|[2"00ll '). ei- 
ther via binary coalescence or collisions. These are the most 
likely to stand out as BSs in the cluster colour-magnitude 
diagram. 

The distribution of binary orbital separations should af- 
fect not only the frequency of mass-transfer events, but also 
the outcomes of dynamical interactions involving binaries. 
It is the closest binaries that are the mo st likely to tidally 
interact and undergo mass-transfer (e.g. iMathieu fc Gellen 
|2009| ). and the probability of a collision occurring during 
1-1-2 and 2-1-2 inte ractions increases w ith decreasing binary 
semi-major axis (jPregeau et al.l |2004| ). Similar arguments 
can also be extended to orbital eccentricity distributions 
that are richer in high-eccentricty orbits. All of this sug- 
gests that the dependences of the various binary parameter 
distributions on total cluster mass could, in principle, play 
a role in driving the observed correlations (or lack thereof). 

4-2-4 The neglect of dynamics 

If a given BS currently resides in the cluster core, this does 
not necessarily mean that it formed there. In particular, 
many BSs could have either formed outside the core be- 
fore migrating in due to dynamical friction, or they could 
have formed inside the core from binary progenitors that 
recently migrated in (e.g. MapcUi ot al. 2006:). We found in 
iLeigh. Sills fc Kniggel (|201ld ) that the observed dependence 
of BS numbers on core mass can only be reproduced if at 
least some BSs did indeed recently migrate in from outside 
the core. Specifically, we found that this contributed to low- 
ering the power-law index on Mcore- This is because the 
time-scale for two-body relaxation increases with increasing 
cluster mass, so that fewer BSs formed outside the cores of 
more massive clusters have had sufficient time to migrate 
in via dynamical friction. If correct, this predicts that the 
global numbers of BSs should correlate more linearly with 
the total cluster mass than we have found for the relation 
between the numbers of BSs in the core and the core masses. 
This can be tested directly using a large sample of cluster 
CMDs derived using a field of view that extends out to the 
tidal radius in all clusters. The on-goin g work of, for exam- 
ple. iFekadu. Sandauist fc Boltd 1I2OO7I). iDalessandro et al.l 



(|2009l), ICarraro fc SeleznevI (|201ll ). iBeccari et all (|201lD . 
and lSanna et al.l (|2012l ) should prove very useful in this re- 
gard. This is because they have slowly been compiling a large 
sample of CMDs with nearly complete spatial coverage, and 
it will be possible to compile from this a homogeneous sam- 
ple of cluster CMDs for which the field of view consistently 
includes the entire cluster. 

On the other hand, most normal binaries currently pop- 
ulating the core are likely to have spent a signi ficant fraction 
of their lives there (e.g. iHeggie fc HutI 120031) . This means 
that they have had plenty of time to have been affected 
by the cluster dynamics. Therefore, any BSs currently in 
the core that recently migrated in from the cluster outskirts 
were formed from a more primordial component of the to- 
tal binary population, whereas the present-day core binary 
fractions reflect a more dynamically-processed component. 
The effects of this could include a weakening of the corre- 
lation between the observed number of BSs in the core and 
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the present-day number of binaries in the core. However, it 
would not affect an underlying correlation with the cluster 
or core mass. This would be consistent with the data, so this 
effect may contribute to our results. 



4-2.5 Selection effecs 

One final effect worth considering is the reliance of our BS se- 
lection criteria on location in the colour-magnitude diagram, 
fn particular, we find that the number of BSs scales sub- 
linearly with the core mass. Therefore, in order for some is- 
sue with our CMD-based selection criteria to explain our re- 
sults, we require that proportionately fewer BSs appear con- 
siderably brighter and bluer than the main-sequence turn- 
off in massive cluster cores when compared to low-mass 
cores. This could arise if, for example, more merger and 
mass-transfer products appear hidden along the MS in pref- 
erentially massive clusters, instead of appearing distinctly 
brighter and bluer than the MSTO. If correct, this predicts 
that the average BS luminosity should decrease with increas- 
ing core mass, ft follows that the average BS mass should 
also decrease with increasing core mass, since previous stud- 
ies have shown that t he luminosities o f BSs are correlated 
with their masses fe.g. lSills et"al]|200lD . This offers a useful 
test of the idea that the number of "blue stragglers" hidden 
along the MS in the CMD depends systematically on the 
total cluster mass. 



5 SUMMARY 

We have carried out a statistical analysis to study the origins 
of blue stragglers in a large sample of Galactic globular clus- 
ters, based on data obtained as part of the ACS Survey for 
Globular Clusters. The main novel ingredient in our analysis 
are empirically estimated core binary fractions, which allow 
us to estimate the number of core binaries, as well as the 
1+2 and 2+2 collision rates. Contrary to our expectations, 
we have found that none of these observationally estimated 
quantities yield correlations with BS numbers that improve 
upon the previously known sub-linear correlation between 
BS numbers and core mass. This is despite the fact that the 
binary fractions themselves anti-correlate strongly with core 
mass, just as expected in a simple binary evolution model, 
where the number of BSs would scale linearly with the num- 
ber of binaries. 

We have explored several possible explanations for our 
results. The simplest, and most appealing, is that observa- 
tional uncertainties affecting the core binary fractions exceed 
the intrinsic scatter of the relationship between binary frac- 
tions and core mass. This could reconcile the data with the 
binary evolution model. In the context of the binary evolu- 
tion model, this would explain why the product of binary 
fraction and core mass is a poorer predictor of BS num- 
bers than core mass alone, and also why the relationship 
between the observed numbers of binaries and BS numbers 
is sub-linear. 

If this explanation is correct, it would imply that core 
binary fractions are tightly coupled to the core or total clus- 
ter mass. This would be of considerable significance for the 
dynamical evolution of globular clusters, and provides an 



important benchmark for simulations attempting to under- 
stand their present-day properties. 
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